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Abstract. Helioseismic inversions for the rotation rate have 
established the presence of a tachocline near the base of the 
solar convection zone. We show that the tachocline produces a 
characteristic oscillatory signature in the splitting coefficients 
of low degree modes, which could be observed on distant stars. 
Using this signature it may be possible to determine the char- 
acteristics of the tachocline using only low degree modes. The 
limitations of this technique in terms of observational uncer- 
tainties are discussed, to assess the possibihty of detecting 
tachoclines on distant stars. 
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1. Introduction 

Helioseismic inversions of observed splitting coefficients have 
enabled u s to s tudy the rotati on ra te inside the Sun (Thomp- 

1998| ). From these inversions it 



X 



son et al. 1996; Schou et al 



has been established that there is a shear layer near the base of 
the convection zone where the rotation rate undergoes a tran- 
sition from differential rotation inside the convection zone to 
almost uniform rotation in the radiative interior. This layer has 
been referred to as the tachocline (Spiegel & Zahn 1992). The 
characteristics of this layer have been studied using helioseis- 
mic data (Kosovichev 199f; Basu 1997, Antia, Basu & Chitre 



1998; Charbonneau et al. 1999; Corbard et al. 1999). Neverthe 



less, the origin of this shear layer is not yet clear and it would 
be instructive to probe the possible existence of these layers 
in distant stars. Such a study would help us in understanding 
the formation of tachoclines and to test the theories of angular 
momentum transport in stellar interiors. 

The solar tachocline has been detected using frequency 
splittings for modes of low and intermediate degree, £. All these 
modes are not expected to be detected on other stars. In order 
to detect a tachocline on distant stars we have to look for the 
signature of a tachocline in low degree modes {£ < 3), which 
are the only modes that can be detected on these stars. It has 
been shown that rapid variations in the sound speed in the stel- 
lar interior, like those arising at the base of the convection zone 



leave a characteristic oscillatory signature in the mean frequen- 
cies of low degree modes (Gough 1990; Monteiro, Christensen- 
Dalsgaard & Thompson 1994; Roxburgh & Vorontsov 1994; 
Basu, Antia & Narasimha 1994). From this oscillatory signa- 
ture it has been possible to put limits on the extent of over- 
shoot below the solar convection zone (Basu 1997). Monteiro, 



Christensen-Dalsgaard & Thompson (2000) have pointed out 
that this oscillatory signature can be used to study the location 
of the base of the convection zone as well as the extent of over- 
shoot below this base in other stars using asteroseismic data for 
only low degree modes. 

In this work we study the signature of a tachocline in the 
low degree modes. Since the tachocline is a narrow layer where 
the rotation rate varies rapidly, we would expect an oscillatory 
signature in the corresponding splitting coefficients. Using a 
model for the solar tachocline, we show that this oscillatory 
signature is indeed present in frequency splitting coefficients 
of the low degree modes. Further, it is, in principle, possible to 
determine the characteristics of the tachocline, like its location, 
width and the extent of variation in the rotation rate across the 
tachocline using this oscillatory signature. 

2. The technique 

The frequency of an eigenmode of a given degree £, radial or- 
der, n, and azimuthal order, ni can be expressed in terms of the 
splitting coefficients, using the expansion 
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(1) 



Here Vni is the mean frequency of the (n, £) multiplet, Cj {n, £) 
are the splitting coefficients and V^^^ jm) a re orthogonal poly- 
nomials in m (RitzwoUer & Lavely 1991). If the rotation en- 
ergy in the star is much smaller than the gravitational energy, 
then rotation can be treated as a small perturbation over the 
non-rotating spherically symmetric star. In this approximation 
the splitting coefficients can be determined by the correspond- 
ing component of the rotation rate. Following Ritzwoller & 
Lavely ( |1991[ ), we express the rotation rate as 
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where is the colatitude, Pj are the Legendre polynomials and 
Wj (r) are expansion coefficients for the rotation rate which are 
related to the splitting coefficients cj {n, i) by 



{n,e)= Wj{r)K'f'^\r)r'^ dr, 



(3) 



where the kernels k''^'''' {r) are defined in terms of eigenfunc- 
tions of the mode (Ritzwoller & Lavely 1991 ). For the Sun it 



is found that a major part of the transition in the tachocline is 
deter mined by the splitting coefficient C3(n, t) (Antia, Basu & 
Chitre 1998 ). The corresponding component of the rotation rate 
is defined by 



1^3 (r,0) 



r sin 9 dO 



P3(cos0) = ^^^^|(5cos2 0-l).(4) 
r 2 



The variation in the other components across the tachocline 
turns out to be an order of magnitude smaller The splitting co- 
efficient C3(n, t) is defined only for £ > 2 and hence only these 
modes can be used to study this component. In the absence of 
any theoretical understanding of formation of a tachocline, it is 
difficult to predict the expected form or magnitude of variation 
for other stars. In this work, we assume that in other stars too 
the dominant variation is in the component. Although, the 
analysis is similar for all splitting coefficients, the higher order 
coefficients can only be determined for modes of higher de- 
gree. Since for other stars we can only expect to detect modes 
with degree I = 0, 1,2,3, it is not possible to study transi- 
tion in higher order terms. On the other hand, if the dominant 
variation in other stars is found in the spherically symmetric 
component corresponding to ci(n, t) it will be easier to detect 
the tachocline, but we will not consider this possibility in this 
work. 

In the asymptotic limit (Christensen-Dalsgaard & 
Berthomieu 



1991 ) we can approximate the kernel as 



c(r) 



where t is the acoustic depth given by 



dr. 



(5) 



From Eq. |] it is clear that if we have a rotation rate which is 
piecewise constant with transition at r = ra, then the corre- 
sponding splitting coefficient will have an oscillatory term of 
the form sin(2ijjr + 0), similar to that observed in the mean 
frequency due to transition in the sound speed (Monteiro et al. 
1994; Basu et al. |1994 ). The oscillatory part arises because the 



integral over the region covering an entire wavelength of oscil- 
lation will give an average value, while the last part between 
the nearest node of the eigenfunction and the location of tran- 
sition, Td will give an oscillatory contribution. Thus we can 
express the splitting coefficient as 



Cj{n,t) =c''\n,t}+A{uj)i 



2lut 



+ 1) 



,(6) 



where c^"'' (n, £) is the smooth part of the coefficient, which 
may arise from possible smooth variations in the rotation rate 
with depth. The amplitude A{uj) is a smooth function of uj, 



while 7 is a constant. The term involving 7 arises from a more 
accurate expression for the vertical wavenumber kr (Monteiro 



et al. 19941. This term can be neglected while considering only 
low degree modes, since it is very small for low £. In general, 
the amplitude also may depend on £, but at low degrees these 
terms are likely to be small. 



Following Basu et al. ( 1994 ), we take the fourth difference 
of the splitting coefficients with respect to n to enhance the 
oscillatory signal. Another advantage of taking the fourth dif- 
ference is that the smooth part of the splitting coefficients be- 
comes negligible and we do not need to include it in our anal- 
ysis. This will of course, depend on the smooth component of 
variation of the rotation rate, but even for a realistic solar rota- 
tion profile this component is found to be negligible. 

In order to illustrate this oscillatory signature we assume a 
model tachocUne rotation profile of the form 



5f}(5cos2 6'- 1) 
1 + exp[(rd - r)/w] ' 



(7) 



where 6il is the extent of variation in the rotation rate across 
the tachocline, w is the half-width of the transition layer and 
Td is its mid-point. It may be not ed tha t this form is differ- 
ent from that used by Kosovichev ( 1996 ) and Charbon neau e t 
al. ( 1999 ), and in particular, as explained by Antia et al. ( 1998 ), 
the definition of width is also different. The half-width w in 
our definition should be multiplied by 4.9 to obtain the width 
as defined by them. Using this model rotation profile we can 
calculate the corresponding splitting coefficients for the Sun. 
Because of the choice of latitudinal dependence only the split- 
ting coefficient 03(71, £) is found to be non-zero and Fig. [l|(a) 
shows the fourth difference of these coefficients. For illustra- 
tion we have included modes with 2 < ^ < 10. The oscillatory 
trend is clearly seen in these differences. The 'wavelength' of 
these oscillations depends on the depth r^, while the amplitude 
is determined by the width and Jfi. Since the splitting coeffi- 
cients and hence the fourth differences scale linearly with Sfl 
we have used only one value 5il — 20 nHz which is the typical 
variation across the solar tachocline. In order to characterize 
the oscillatory signal we fit the following form to these differ- 
ences: 



d C3{n,£) = [ao -\ ! 7 sm 2ti)T + 



(8) 



The parameters ao, ai, 02, r, 0, 7 can be determined by a non- 
linear least squares fit. 

For distant stars it will be possible to detect modes with £ = 
0,1,2,3 only and hence the number of modes will be highly 
restricted. Even then it is possible to fit the oscillatory part and 
Fig.|(b) shows one such fit. In this case the parameter 7 is not 
relevant as the corresponding term is very small and we fit only 
the five parameters ao, ai, 02, r and (p. All the results presented 
in this paper are obtained using modes for £ = 2,3 only, unless 
mentioned otherwise. 
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Fig. 1. Fourth difference of the splitting coefficient C3, for 50, — 20 nHz, w — O.OO5i?0 and ra = O.Ti?©. The continuous line 
shows a fit according to Eq. (a) Modes with 2 < £ < 10 have been included. The frequency has been shifted to account for the 
^-dependence of the term in sine function, (b) Only modes with £ = 2 and £ — i have been included. 



3. Results 

In order to study the signature of the tachocline in the splitting 
coefficients we calculate c^{n, £) with rotation profiles defined 
by Eq. ^ using different values of half- width w, and position 
of the tachocline. We fit the oscillatory form defined by Eq. |] 
to each of these data sets. It turns out that the position of the 
tachocline affects mainly the parameter r, which is close to the 
acoustic depth of the transition layer. Therefore, in this work 
we have only shown results with = O.7i?0. Conversely, the 
location of the tachocline may be determined from the parame- 
ter T. In order to study the effect of width we use different val- 
ues of ui and the fitted amplitude (j4(cj) = ao + ai/uj + a2/uj'^) 
is shown in Fig. ^ It can be seen that in all cases the amplitude 
decreases with increasing uj, which may be expected as the 
modes with higher frequency have larger number of nodes in 
radial direction and hence the radial wavelength will be smaller, 
thus giving a smaller contribution to the integral in Eq. |[ At 
larger width the amplitude reduces and the tachocline profile 
also gives a contribution to the smooth part at low frequencies 
and it is not possible to use the simple form to fit the data. It 
turns out that the extent of variation in the amplitude across 
the frequency range depends on the width of the tachocline. 
For small width the variation is smaller, while with increas- 
ing width, the variation in amplitude increases. For example, 
for w = O.OO2i?0, 0.003i?Q, 0.0057?© and O.Oli?©, the ra- 
tios of amplitudes at the two limits in Fig. ^ are 2.5, 3.0, 5.5 
and 100 respectively. Thus, if we have data covering the entire 
frequency range it should be possible to determine the width 
using the extent of variation in amplitude. Once the width is 
determined, we can determine Sfl from the amplitude of the 
oscillatory term at the low frequency end. Thus we should be 
able to determine the characteristics like position, width and 
Sn of the tachocline using only low degree modes. 

In all the foregoing calculations we have used the exact 
sphtting coefficients as calculated for a given rotation profile. 
In real data, there would naturally be errors associated with 
each splitting coefficient. In order to simulate real data we have 
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Fig. 2. Comparison of amplitude of oscillatory signal as a func- 
tion of frequency for different widths of the transition region 
for Sn = 20 nHz and = 0.7Rq. 



constructed artificial data sets where random errors with stan- 
dard deviation a are added to all the splitting coefficients. For 
simplicity, we assume that error is same in all these coeffi- 
cients. Using 100 different realizations of errors we can esti- 
mate the expected errors in each of the fitted parameters, the 
results of which are summarized in Table |l]. For low errors 
(cr ^ 0.02 nHz) it is indeed possible to determine all the param- 
eters to reasonable accuracy and the eiTor in each parameter is 
proportional to the assumed error in the splitting coefficients. 

However, as error increases it is not possible to determine 
the amplitude variation correctly, but if we fix ai =0 and 
02 = then we can still determine the mean amplitude of 
the oscillatory part. For Sfl — 20 nHz, we find that it is pos- 
sible to detect the oscillatory signal reliably when the error 
(7 ^5 0.2 nHz, which corresponds to an error of about 1.6 nHz 
in the fourth differences. If we assume that the tachocline in 
other stars would also be located close to the base of the con- 
vection zone, then its depth can be determined by lookin g at the 
oscillatory part in the mean frequencies (Monteiro et al. 20Q(\ ) 
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Table 1. Values of fitted parameters in oscillatory signal for 



= 20 nHz, w 


= O.OO5i?0, 


ra = O.7i?0 




Error (a) 


ao 


T 




(nHz) 


(nHz) 


(sec) 


(rad) 


0.010 


0.68 ±0.02 


2323 ± 5 


4.70 ±0.14 


0.020 


0.68 ± 0.04 


2324 ± 10 


4.68 ±0.28 


0.050 


0.70 ±0.10 


2331 ± 32 


4.59 ±0.81 


0.100 


0.81 ±0.19 


2330 ± 135 


4.53 ± 1.40 


0.200 


1.25 ± 0.45 


2290 ± 242 


4.79 ± 1.63 


0.200 


0.79 ± 0.36 


2323 


4.70 ± 0.78 



which determines the base of the convection zone. In that case 
we can fix r to the value determined from the convection zone 
base to improve the detectability of the oscillatory signal. The 
last line in Table [l| shows the results of simulations obtained 
when r was kept fixed. It is clear that the estimated errors in 
this case are smaller. If we have the splitting coefficients for 
only £ = 2 modes then the errors are found to be \/2 times 
larger For the solar case we can use higher degree modes also, 
but even in that case the error in the observed splitting coeffi- 
cients is slightly larger than this limit and it is not possible to 
detect the signal reliably. 



4. Discussion 

We have shown that a sharp transition in the rotation rate that 
is expected in the tachocline region gives rise to an oscillatory 
signal in the splitting coefficients of the low degree modes. The 
'wavelength' of the oscillatory signal is determined by the po- 
sition of the tachocline, while the variation in amplitude across 
the frequency range is determined by its width. The amplitude 
of signal is of course proportional to the extent of the variation 
in the rotation rate. Thus from this oscillatory signal it is, in 
principle, possible to determine the position, width and Sft for 
the tachocline. We are looking for the signal in odd splitting 
coefficients, which can arise only due to rotation, and not from 
magnetic fields or structural variations in the stellar interior 
In this discussion we have not included any possible latitudi- 
nal variation in the characteristics of the tachocline. The solar 
tachocline is known to be prolate (Charbonneau et al 



1999) 



and this variation would also affect the oscillatory signal in low 
degree modes. Using only low degree modes it may not be pos- 
sible to disentangle all variations in the tachocline, but so long 
as the latitudinal variation is small, as is the case for the Sun, 
the mean properties of the tachocline can be determined from 
the low degree modes. 

Since the amplitude of the oscillatory signal is very small, it 
will be necessary to find the splitting coefficients accurately to 
determine the characteristics of the tachocline. From our sim- 
ulations it appears that the required accuracy in the splitting 
coefficients is ^ 0.2 nHz for Sfl = 20 nHz. Even for the Sun, 
this level of accuracy has not been achieved so far and we do 
not expect it to be achieved for other stars in near future. But 
there is no reason to assume that the variation in rotation rate 
in all stars will be comparable to that in the Sun. In particu- 



lar, for stars which are fast rotators, we would expect Sfl to be 
correspondingly larger Even for stars with similar rotation rate 
there may be some variation in or in the amplitude of oscil- 
latory signal for the same 6il. In this work, all the calculations 
have been done for a solar model. In other stars the amplitude 
may be somewhat different. If the star is rotating very rapidly, 
the linear approximation used to relate the splitting coefficients 
to the rotation rate may not be admissible. But we would ex- 
pect that for a star which is rotating about 50 times faster than 
the Sun, this approximation may still be applicable and in that 
case an accuracy of about 10 nHz may be sufficient to detect 
the oscillatory signal. For stars with larger M/B?, this limit 
may also be larger. Similarly, if we choose stars with larger 
differential rotation, or with favourable amplitude, this number 
may go up further. Even if stars are rotating more rapidly the 
oscillatory signal may still be present, but mode identification 
and interpretation may be more difficult. The upcoming aster- 
oseismic missions like COROT (Baglin et al. | 1998|), M OST 



(Matthews |1998| ) and MONS (Kjeldsen &. Bedding |199§) have 
a planned frequency resolution of 100 nHz, which may not be 
sufficient to detect this oscillatory signal. But with some im- 
provements in instruments and longer observations of a few se- 
lected stars, which are known to be rotating fast and preferably 
having larger differential rotation, it may be possible to detect 
this signal in not too distant future. 
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